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Abstract 

We study a semilinear PDE generalizing the Fujita equation whose evo- 
lution operator is the sum of a fractional power of the Laplacian and a con- 
vex non-linearity. Using the Feynman-Kac representation we prove criteria for 
asymptotic extinction versus finite time blow up of positive solutions based on 
comparison with global solutions. For a critical power non-linearity we obtain 
a two-parameter family of radially symmetric stationary solutions. 

By extending the method of moving planes to fractional powers of the Lapla- 
cian we prove that all positive steady states of the corresponding equation in a 
finite ball are radially symmetric. 

Resume 

Nous etudions une equation de reaction-diffusion semilineaire (generalisant I'equa- 
tion de Fujita), dont I'operateur d'evolution est la somme d'une puissance frac- 
tionnelle du Laplacien et d'une non-linearite convexe. A I'aide de la representation 
de Feynman-Kac nous exhibons des criteres entrainant I'extinction asympto- 
tique, respectivement I'explosion en temps fini, de solutions positives. Ces 
criteres s'obtiennent en comparant avec des solutions globales. Pour une cer- 
taine puissance critique de la non-linearite nous obtenons une famille parametrisee 
de solutions stationnaires a symetrie radiale. 

Par extension de la methode de deplacement d'hyperplans a des puissances 
fractionnelles du Laplacien, nous prouvons que toute solution positive station- 
naire de I'equation correspondante dans une boule finie comporte une symetrie 
radiale. 
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1 Introduction 



We consider the "generalized Fujita equation' 



d_ 

di 



m(0, x) 



Lu{t, x) + G{u{t, x)), t > 0, X e R"* 
> 0, 



(1.1) 



where L = Aq, is the fractional power — (— A)"/^ of the Laplacian, < a < 2, 
and G : R+ — R+ is a convex function satisfying conditions (2.1) and (2.2) below. 
Solutions will be understood in the mild sense (see e.g. [20]) so that (1.1) makes sense 
for any non-negative bounded measurable function ip on R"'. 

A well-known fact is that for any non-trivial initial value (p there exists a number 
G (0, oo] such that (1.1) has a unique solution u on R'^ x [0, T^) which is bounded 
on R'^ X [0,T] for any < T < T^, and if < oo, then ||m(-, t)||Loo(Rd) — >■ oo as 



When = oo we say that u is a global solution, and when T<^ < oo we say that 
u blows up in finite time or that u is non-global. 

The study of blow up properties of (1.1) goes back to the fundamental work of 
Fujita [11], who studied Eq. (1.1) with a = 2 and G{z) = z^+f^, /3 > 0. The 
investigation of (1.1) with a general a was initiated by Sugitani [22], who showed 
that ii d < a/ (3, then for any non- vanishing initial condition the solution blows up 
in finite time. Using a Feynman-Kac representation for the solutions of semilinear 
problems of the form (1.1), this conclusion was re-derived in [4] and the corresponding 
behaviors of equations with time-dependent nonlinearities and of various systems of 
semilinear pde's were studied. 

It is known (e.g. [18], [19]) that in supercritical dimensions d > a/P, Eq. (1.1) 
admits global as well as non-global positive solutions, depending on the "size" of 
the initial condition. For short we address this parameter constellation as the global 
regime. 

In the first part of the present note we prove two comparison criteria in the global 
regime: 

(i) Assume the initial value (p > leads to a globally bounded solution. Then any 
initial value ip with < < {1 — e)ip, e > 0, gives rise to a solution converging to 
zero. 

(ii) Assume the initial value 9? > leads to a solution which is uniformly bounded 
away from for alH > and all x in some ball C R'^. Then any initial value if) with 
ip > {1 + e)ip, e > 0, gives rise to a solution which blows up in finite time. 

The essential tool in proving (i) and (ii) is the probabilistic representation of the 
solution of (1.1) provided by the Feynman-Kac formula, that was obtained in [4] (see 
(2.4) below). 

Natural candidates for the comparison in (i) and (ii) are (time-) stationary solu- 
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tions of (1.1), i.e. solutions of the "elliptic" equation 

A„n(x) + G{u{x)) = 0, X e R'^. (1.2) 
In the case a = 2 and G{z) = ^^+^, it is known that (see [14], [9], [12], [23]) 

- for d > 2, 1 + /5 < {d + 2)/{d — 2), apart from u = 0, no bounded non-negative 
solution of (1.2) exists 

- for d > 2, 13 = {d + 2)/{d — 2) — 1, all bounded stationary solutions of (1.2) are 
given by the family 

UcA^) = - A{d{d - 2)Y''-'y' cxeR^ AeR+. (1.3) 

(did - 2) + {A^/(d~^)\\x - c\\f] 

(Note that the two parameters of the family are the symmetry center c of u and its 
value A at c.) 

- ioT d > 2, P > {d + 2)/{d — 2) — l, there exists a one-parameter family ua, A > 
of solutions of (1.2) with the properties: ua{x) is symmetric around x = 0, 

ua{0) = A, llxf/f^UAillxW) K{d,(3) as ^ oo 

where K{d, (3) is a constant not depending on A. 

In the case a < 2, much less is known. For d > a, (3 = {d + a) / {d — a) — 1, 
we specify in Proposition 3.1 a two-parameter family Mc,a, c G R*^, < A < oo, of 
radially symmetric solutions of 

A^u{x) + u^+^{x) = 0, xeR'^ (1.4) 

with the property 

'Wc,a(c) = A, ||x||'^""Mc,A(||a:||) — ^ K(d,a,P) as — > oo 

where K{d,a,l3) is a positive constant. A natural conjecture now is that, like in 
the Laplacian case, for /? = {d + a)/ {d — a) — 1 the Uc,a constitute all the bounded 
solutions of (1.4), and that for (3 < {d + a)/{d — a) — l there are no bounded non-zero 
solutions at all. 

As a first step to answer these questions, in Section 4 we make use of the so 
called "method of moving planes", which is well known in the Laplacian case, to 
show symmetry of positive solutions in a ball. Essential tools (like Hopf's boundary 
lemma) can be carried over to the a-Laplacian case. With this method, we were able 
to show (see Theorem 4.1) that the following equation has only radially symmetric 
solutions: 

Aau{x) + F{u{x)) = 0, xeB, 
u{x) = 0, xeR'^\B. 

Here B is an open ball centered around and F : R^ R_|_ is non- decreasing. We 
conjecture that an analogous statement is valid for F[z) = with p = [d+a) /{d — a) 
and R*^ instead of B. In Section 4.3 we describe several problems we think one would 
have to overcome for carrying over the moving planes method to the unbounded space 
setting in this case. 
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2 Two comparison criteria 

In this section we assume that the function G in Eq. (1.1) satisfies the conditions 



and 

dz 



G{z] 



< oo (2.2) 



for certain positive numbers (3 and 6. 



Lemma 2.1 Let G be a convex function satisfying (2.1), and £ > 0. For any M > 

there exists 5' > such that 

^((l+e)^) , „Giz) , 

> 1 + e')^^ forO<z< M. 

{l + e)z ^ ' z 



Proof By considering G/c instead of G we can assume that c = 1. Given e > 
there exists 5 > such that 

[l-e)z^ <^<{l + e)z^ 
for z e (0, 5). Take €<{{! + - 1)/((1 + + 1). Then for z < 5/(1 + e) := xq, 
''((I + > (1 - .1(1 + efrf > (IziHlif)!^ = (1 + 



{l + e)z ' ' (1 + e) 

where c' > 0. Since z ^ G{z)/z is continuous and strictly increasing in (0, cxo) it 
follows that inf,6[^„,M] /^) > 1 + c" with c" > 0. Taking = c' Ac" yields 

the assertion. □ 

Let us observe that if w is a globally bounded solution of (1.1) we necessarily have, 
for all X, 

Ptv{0,x)^0 as t — > cxD, 
where (Pt) is the semigroup with generator L. Indeed, from the integral form of (1.1) 

v{t,x) = Ptv{0,x)+ / Pt^sG{v{s, x)) ds 

Jo 

> Ptv{0,x)+ I Pt^sGiP,v{0,x))ds 



> Ptv{0,x)+ [ G{Ptv{0,x))ds 

Jo 

> tG{Ptv{0,x)), 



4 



where we used in the first inequahty that Ptv{0, ■) < v{t, ■), and Jensen's inequahty 
after the second line. It follows from the global boundedness of v that 

lim Ptv{0, x) < lim G"^(Const. r^) = 0. (2.3) 

i— »oo t— >oo 

Proposition 2.1 Let G be a convex, increasing function satisfying (2.1) and (2.2). 
Assume the initial value > leads to a globally bounded solution of (1.1). Then 
any initial value ip with < < {1 — e)ip, e > 0, gives rise to a solution converging 
uniformly to zero. 



Proof Recall that the Feynman-Kac representation of solutions of (1.1) is given by 
(see [4]) 



u{t, x) 



uiO,y)ptiy,x)Ey 



exp 



G{u{s,Xs)) 



ds 



Xt = x 



dy, 



(2.4) 



where (Xt) is the Levy process with generator L, and pt{x,y), t > 0, x,y E R*^, are 
its transition densities. 

Suppose that v is a globally bounded solution of (1.1) and that < u(0, ■) < 
(1 — e)v{0, ■) where < e < 1. As (1.1) preserves ordering we have u{t,x) < v{t,x) 
for alH > and x G R"^, which together with (2.4) improves to 



u{t,x) < / {1- e)v{0,y)pt{y,x)Ey 
= (1 — e)v{t, x) 



exp 



G{v{s,X,)) 
v{s,Xs) 



ds 



Xt = X 



dy 



uniformly in t and x. Inserting this bound again into the Feynman-Kac representation 
of u yields 



u{t,x) < / u{0,y)pt{y,x)Ey 



exp 



G{{l-e)v{s,Xs)) 
(l-e)v(s.X^ 



ds 



Xt = X 



dy. 



Putting 2;(t, x) := (1 — e)v{t, x) in the above inequality and using Lemma 2.1 (with e 
in Lemma 2.1 substituted hy e := e/ [1 — e)) we get 



u{t,x) < / u{Q,y)pt{y,x)Ey 



exp 



* 1 G{{l+e)z{s,Xs)) 



l + e' {l + e)z{s,X,) 



ds 



X, = x 



dy 



uiO,y)ptiy,x)Ey 



R'' 



Xt = X 



dy, 



where e' > is given by Lemma 2.1, e" := e' /{I + e') and 

G{v{s,X,)) 



A 



t ■- 



v{s,X,) 



ds. 



(2.5) 
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Thus, 



Xt = x 



dy 



+ I u 



Since e^* > (Pit;(0, x))-^/^ implies e^^"^")^* < e^* (Pt^;(0, x))^"/^ we obtain 
M(t,x) < / m(0,?/)pj(?/,x) (Ptt;(0,x))"^^^ 

+ (Pii;(0, x))^"/' /" n(0, y)vt{y. x)E^ [e^* | = x] dy 



< il-s) (^iPMo,x)Y/' + v{t,x)iPMo,x)r"/' 

which tends to uniformly as t ^ cxd due to (2.3). 



□ 



Proposition 2.2 Let G be a convex, increasing function satisfying (2.1) and (2.2). 
Assume the initial value = w(0, •) > leads to a globally bounded solution of (1.1) 
which for some open ball B C R'^ and some k > obeys 



inf v{t^x) > K, for all sufficiently large t > 0. 



(2.6) 



Then for any e > 0, the initial condition (1 + e)(p leads to blow-up in finite time. 
Proof By the Feynman-Kac formula, 

v{t,x) = / v{0,y)pt{y,x)Ey [e^*| Xt = x] dy, 
where At is given by (2.5). If K > then 

/ viO,y)ptiy,x)Ey [e'^'-At < K\ Xt = x] dy < e^E,t;(0,Xt) ^ 



as t ^ oo uniformly in x due to (2.3). Therefore, for all > there exists Tq 
Tq{K, 7) > such that for t > Tq 

v(t, x) 



< 2. 



(2.7) 



J^,v{0,y)ptiy,x)Ey[e^^-At>K\ Xt = x] dy 

Without loss of generality we can assume that inf^.gBi(o) ^('f? 2;) > k, for all t large 
enough, where Br{x) denotes the ball in R'^ of radius r centered at x. Arguing as 
above we check via the Feynman-Kac representation that for all t > and x G R*^, 
u{t,x) > (1 + e)v{t,x). Plugging this again into the Feynman-Kac representation for 
u yields 

Giil + e)vis,X,)) 



u{t,x) > {l+e) j v{0,y)ptiy,x)Ey exp / ^j- 
jRd L Jo U 



+ e)v(s,Xs 



ds 



Xt = X 



dy 



> (l+e) v{0,y)pt{y,x)Ey 

jRd 



Xt = X 



dy 



for some e' > by Lemma 2.1. Using this and (2.7) we obtain for given K > and t 
big enough that 

u{t,x) > (l+£)e^"'/" v{0,y)pt{y,x)Ey[e^';At>K\Xt = x]dy 

Hence, for any > we find inix^Biio) u{t, x) > K{l+e)e^^' /2 for all sufficiently large 
t. As is well known (see e.g. [16]), this inequality together with (2.2) are sufficient 
for finite-time blowup of u. □ 

Corollary 2.1 Let G be a convex, increasing function satisfying (2.1) and (2.2), and 
ip > a non-trivial positive bounded solution of 

Lip{x) + G{ip{x)) = 0. 

a) For each e > the solution u of (1.1) with initial value u{0,x) = (1 + e)ip{x), 
X G R'^, blows up in finite time. 

b) For each s G (0, 1) the solution u of (1.1) with initial value m(0, x) = {l—e)(f{x), 
X G R'^ converges uniformly to as t oo. 

Proof This is immediate from propositions 2.1 and 2.2. □ 



3 A class of radially symmetric stationary solu- 
tions 

We now set out to specify a family of positive stationary solutions of (1.1) in the 
particular case of G{z) = z^, where p = {d + a)/{d — a). Before doing this, we still 
consider the case of a general p, and note that the "elhptic" equation 

Aau{x) + uP{x) = 0, X G R'^ (3.1) 

can be rewritten in integral form as 



u{x) 



POO 

/ E,. [uP{Xt)] dt, X G R^ (3.2) 

^0 



where (Xt) denotes the (symmetric) a-stable process in R'^. Hence, for d < a, due 
to recurrence of (Xt), the only non-negative solutions of (3.1) are u = and u = oo. 
Therefore, we henceforth assume that d > a, in which case (3.2) rewrites as (see [5], 
p. 264) 

''W^/-l^1F^*• '''''' 

Jw \\y ^ "''II 

where A{d, a) := r(i(c/ - a))/[r(ia)2"7r'^/2]. 
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Proposition 3.1 If p = {d + a) /{d — a) , then for any A G (0, oo) and c G R"' the 
function 



solves (3.1). 



A 



1 + (^^2/(d-a)2-i (r [d^) /r (^) )"'^" ||x - C||)' 



{d~a)/2 



Proof Without loss of generality we assume that c is the origin. Due to (3.3) it 
suffices to show that 



Uo,a{x) 



A{d,a)uljs,{y) 



I d—a 



dy, X G 



\\y-x\ 

Let us write a := A-2/('^-")2 (r (^) /V (^) )^^". We ffist note that 

A 



(3.4) 



Uo,a[^) 



1 +47r2||-^||2)(d-a)/2 
II 2-7™, II ' 



AB, 



2'Ka 



(see [10], p. 155), where for any / G L-'^(R"), f{x) := j-^^ e ^'^^^'^/(y) is the Fourier 
transform of /, and for any complex w with Re(-u;) > 



Vm^nY'^ Jo 



.{w-d)/2-l^-T-\\x\\^/iT 



Hence 

uoa{x) = A 
Notice that 



B, 



d—a 



2TXa 



(x) = A{2-KaYBd-c{'2Tiax) = A{27raYBd-ai-27rax). 



Bw{x) 



2{d-w)/2+l 

r (f ) (47r)'^/2 



\x\ 



^{w~d)/2^ 



{d-w)/2[.\\x\\), 



xen" 



(3.5) 



where for any complex z/, is the Macdonald's function ([24], §6-22), also known as 
modified Bessel function of the second kind, which is given by 



= \ (i. 



It follows that 



°/2K«/2fll27rax| 



r(^) 



(3.6) 



To compute the Fourier transform of the other side of (3.4) we use the convolution 
theorem, (3.5) and 



[A{d, 



a 



~{d-a) 



]{x) = (27r||x||)~° , X G R*^, < Re(a) < d, 
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e.g. [10], p. 154) to obtain 



Rd \\y 



I d—a 



■dy 



{x) = (27r)-"||x||~"AP 



{d+a)/2 



(271) 



B. 



d+a 



2na 



ix) 



{2iTy\\x\\^''AP{2iTaYBd+a{-27iax) 
2-V'^""^/'||x|r"/'AV+"/' _ „ . K^^/2{\\2nax\ 



r(^) 

Since K^, = K_y ([1], Formula 9.6.6), by comparing the RHS of the last equality with 
that of (3.6) we see that they are equal for the value of a stated at the beginning of 
the proof. The result follows from uniqueness of Fourier transforms. □ 

Remarks 1. Recall [8, 15] that for < a < 2 and d > a the Kelvin transform of u 
is defined by 

(3.7) 



v[x] 



\x 



A simple calculation shows that for any fixed c G R"' the family of solutions {uc,a\a>o 
rendered by Proposition 3.1 is invariant under the Kelvin transform with center c 



Vc{x) :-- 



X — c 



I d—a 



X — C 



\x — c 



2. Moreover, if u is any given regular positive solution of (3.1) (where < a < 2, 
p > 1 and d> a) and x 7^ 0, then its Kelvin transform v{x) satisfies 



^ryVix) 



v'^{x) 



0. 



(3.8) 



||2;||(rf+a)"P{d-a) 

Indeed, let Ga denote the Green's operator corresponding to Aq,, and x 7^ 0. Then 
vP{x) \ f A{d,a) 1 



— Ga 



\rf,\\ {d+a)—p{d—a) 



^j, a ||y|| (d+a)— p(d— o) 



1 



I {d—a)p 



y 



\y\ 



dy 



A{d, a) 
A{d, a) 
A{d, a) 



\\4 



I d~a 



■ u^(z) dz 



X ■ \\z\ 



d~ct 



■ u^z) dz 



u^{z) dz 



\X 



I d—a 



U 



X 



\X\ 



\x\\ ■ z 



v{x) 



d—a 



|a;|| ■ z\\ in the fourth 



where we used the elementary identity — ||fj||| = ||||fjj 

equality. 

3. Proceeding as in the proof of Proposition 3.1 one can verify that a singular 
explicit solution to (3.1) is given by 



2" (r(^) /r(^))^ 



||2;||(rf-a)/2 ' 

and that fixed point of the Kelvin transform (3.7). 



4 Rotational symmetry of solutions in a ball 

Let F : R_|_ R+ be non-decreasing, not identically constant, and a G (0,2). Let 
u : R'^ — s> R_|_ be a non-negative bounded solution to 

AaU + F{u) = 0, X e Bi{0), u = OmBi{OY, (4.1) 

i.e., for all x in the unit ball we have 

u{x) = [ [FiXt); sup,<, < l] dt = [ G^x, y)F{u{y)) dy, (4.2) 

Jo JBiiO) 

where X is the symmetric a-stable process and Ga{x, y) is the corresponding Green's 
function for the unit ball. 

That (4.1) has non-trivial positive solutions follows from [7] at least in the case 
a = 2. Without loss of generality we assume that u ^ 0. Our aim here is to show 
that 



Theorem 4.1 u is rotationally symmetric about the origin. 



Our approach is based on the method of moving planes, a device that goes back to 
Alexandrov [2] and has by now a venerable history in the study of symmetries of 
solutions of pde's, see also [21] and [12]. The idea is as follows: Choose any direction 
in R'^, wlog the Xi-direction, and show that u is mirror symmetric with respect to the 
hyperplane through the origin with this given direction as a normal vector. In order 
to achieve this let us define for A G ( — 1,1) 

Ta := {x G R"^ : Xi = A}, Sa := {x G R'^ : Xi < A} 

and for x G R^ let x'^ := (2A — Xi, X2, . . . , Xd) be the image under reflection along Ta; 
see Figure 1. 

Define the set A by 

A:= |ag (-1,0) :m(x^) > m(x) Vx G Sa, t^u(x) > Vx G Ta H 5i(0) 1 . (4.3) 
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Figure 1 



Observe that by the minimum principle we have u(x'^) > u{x) for x G Sa ni?i(0) and 
A G A, i.e., A G A means that reflection along T\ (strictly) increases the value of u. 
In Section 4.2 we prove that 

supA = (4.4) 

so that by continuity u{—Xi,X2..-,Xd) > u{xi,X2 ■ ■ ■ ,Xct) whenever Xi < 0. By 
considering A > and working in the opposite direction we can then conclude the 
reversed inequality and hence obtain the desired symmetry. 



4.1 Some preparatory lemmas 

Lemma 4.1 A bounded solution u of (4.2) satisfies u G C(R'^) fl C°°(i?i(0)). 



Proof Using the explicit form of the Green's kernel of the ball (see e.g. [8], formula 
(2.3)) one easily checks that interchange of integration and differentiation is justified. 

Lemma 4.2 (Minimum principle) Let D C R*^ he a hounded domain. Suppose 
u : R'' R+ is continuous, AaU < on D, and satisfies u = on D'^. Then either 
u = oru>0 on D. 



Proof Let D, := {x e D 



u[x) 



> e}. By continuity is open. Assume that 



u{Xo) - /q^^ Aau{Xs) ds is a P^ 
= inf{s : Xg G D^}. For each x G -D 



% for some e > 0. Let {Xt) be the a-stable process, and r := inf{s : Xg ^ 
the hitting time of D'^. Then Mt := M(XtAr) - 
martingale for each x G -D. Let furthermore r' := 
we have E^-M^/ 



0, or 



u{x) = [u{X^'^r)]+'E^ 



t'At 



-A,m)(X, 



ds 



> E, [uiXr 



'At) 



> eP^ir' < r) 



> 
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Figure 2 



because (Xt) hits any open subset of D with positive probabihty before exiting from 
D. □ 

We will have occasion to consider the behavior of u at the boundary of the ball. 
In this respect, the following lemma is helpful: 

Lemma 4.3 (Hopf's a-stable boundary lemma) Let D C R*^ be open, u : R'^ 
R+ continuous with u = on D^, A^m < on D, u not identically zero. Let 
Xq G dD satisfy an interior sphere condition, i.e. there exists a ball Bs{xi) C D with 
Bs{xi) n = {xq}, and let v be an outward pointing unit vector at xq- Then 

— m(xo) < 
(in fact, lim£\o(w(xo) — u{xq — eu))/e = —oo). 

Proof Because of the interior sphere condition at Xq we can find a ball Bs{xi) C D 
such that Bs{xi) D D'^ = {xq} and also B C D, where B is the "left half" of a 
spherical shell around xi with interior radius 6 and exterior radius 6' > 6; see Figure 
2. Observe that m > in D by Lemma 4.2, in particular inf^w > because B is 
compact and u continuous. Let r := inf{i : Xt ^ Bs{xi)}. Then 

>E,u{Xr) 

for X G Bs{xi). Take x = Xq — eu with e > small enough, and denote by 7 the angle 
between xi xq and xxq- Then 7 G (— 7r/2,7r/2) because u is an outward pointing 
vector, hence C0S7 > 0. The cosine theorem gives \x — xip = 5^ + — 2(5£:cos7. 
Using the explicit form of the Poisson kernel for the complement of a ball (see e.g. 
[8], formula (2.2) or [6], Thm. A and rescale) we can estimate 

u{x) > E^u{Xr)> IP{x,y)u{y)dy 

Jb 



u{x) = E^ 



u{Xr 



{Ao,u){Xt)dt 
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B J Jb 



y — Xi\'^ — 6 



> c{6'-\x-xr\Y^'>C'e'^/'. 
Thus we see that hmsup£^Q(u(xo) — u{xo — eu))/e = — oo. □ 

Lemma 4.4 Let w : R*^ ^ H be continuous, bounded with w > on Sq, w{x^) = 
—w{x). Let x^ G To be such that there exists a 6 > with AaW < on Bs{x^) fl Sq. 
Then either 

d 

w = 0, or w > on Bs{x^) fl So and — — w{x^) < 0. 

OXi 

Proof Assume w ^ 0. Let x G Bs{x^) fl Sq. If w{x) = we would have (see e.g. [8]) 

\yr 

by the non-triviahty and symmetry of w, in contradiction to the assumption. 

To show that the derivative is non-zero choose 6' > such that sup2Q\B^,(a;,) w > 0. 
Let (Xt) be the a-stable process, r := inf{t : Xt ^ Bs'{x^:)}. Define v{x) := Fj^w^Xr), 
v{x) := Ex Jq {—Aaw){Xs) ds. Observe that 

AaV = in and v = w in R"' \ 

AaV = AaW in Bs'{x*) and = in R"' \ Bs'^x^,). 

Uniqueness of the Dirichlet problem for in Bs'{x^) thus gives w = v + v. We have 
v{x) = J\y_x P{x, y)w{y) dy for x G Bs'{x^:), where the Poisson kernel is given by 
(see e.g. [8], formula (2.2)) 



P{x,y) = a 



5'"^ — \x — 
\y — — 



a/2 



\x-y\ X e Bs'{x^),y ^ Bs'{x^). 



One checks that -^P{x^,y) < for y G Sq and -^P{x^,y) > for y G —So. The 
interchange of integration and differentiation is justified because w is bounded, so we 
can compute 

d f d 

—v{x^)=l —P{x^,y)w{y)dy<0. 

O^l J\y-x^\>5' OXi 

Furthermore, for x G Bs'{x^) 



v{x) = G(x,y){-Aaw){y)dy, 

where the Green kernel for i?5/(x*) is given by (see e.g. [8], formula (2.3) or [6] and 
consider the obvious scaling properties of a-stable processes) 

/Wsi{x,y) a/2-1 
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where ws'{x, y) = (5'^ — |x — a;*p)((5'^ ^ \y ^ x^,\'^)/\x — yp. Inspection shows that for 
X, y G n So we have G{x,y) > G{x,y^), hence > in Bs'{x^,) fl So and by 

symmetry -u < in ^^/(x*) fl (—So). We conclude that {d/dxi)v{x^) < and thus 

d d 
-w{x^) < - — v{x^) < 0. 



dxi dxi 

□ 

4.2 Proof of (4.4) 

We proceed in three steps and show that 

1. A D (—1, —1 + e) for some small e > 0, 

2. A is open. In particular for A G A there exists e > such that [A, A + e) C A. 

3. From 1. and 2. we conclude that A = (— l,Aniax). We finally show that 

Amax 0' 

Step 1. Obviously (—1, 0, . . . , 0) is an outward pointing direction for each x G 9i?i(0)n 
S_i/2. By the boundary lemma and the fact that {d/dxi)u is continuous in -Bi(O) 
there is an open neighborhood D of dBi{0) fl S_i/2 such that {d/dxi)u > on 
D n 5i(0). Choose e > so small that S_i+e n 5i(0) C D. Then -1 + e/2 G A. 

Step 2. We argue by contradiction. Assume there was A^, G A and also a sequence 
(A„) C (—1,0) \ A with A„ \ A^,. iFiom. the definition of A, possibly passing to a 
suitable subsequence (which we again would denote by (A.„)) we can always arrive at 
one of the following possibilities: 



a) There exists a sequence (x„) C -Bi(O), x„ G S;^^, with x„ — > x^, G -Bi(O) and 
u{xn) > u{x^") for all n, or 



b) There exists a sequence (x„) C -Bi(O), x„ G Ta„, with x^ — > x* G -Bi(O) and 
{d/dxi)u{xn) < for all n. 

Assume a) was true. We cannot have x,,. G Sa, because u is continuous and m(x^*) > 
u{x) for X G Sa^ by the above remark. Hence x* G Ta^ fl -Bi(O). But then we have 
{d/dxi)u{x^) = lim„_^oo(M(a;^") - u(x„))/(2d(x„, Ta„)) < 0. By Hopf's boundary 
lemma, this forces x* to be away from dBi{0), but then we obtain a contradiction to 
A^ G A. 



If b) was true we would again find a point x* G Ta, n5i(0) with {d/dxi)u{x^,) < 
and arrive at a contradiction. 

Step 3. iFrom the preceding steps we know that A = (—1, Amax). If m(x'*'™="') = u{x) 
for X G Sa^^^ then we have found a symmetry center. As u is continuous, m = on 
dBi{0) and strictly positive inside -Bi(O) this can only be true for A^ax = 0. Indeed, 
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if Amax < then by continuity we would have u{x) < for x G S^^^^, but with 

u{x) ^ u{x^'^'"'). Define w{x) := u{x^'^'''') — u{x). Observe that w is continuous and 
bounded, non-negative in S^^^^ and w{x^'^''^) = —w{x). For x G S^^^^ fl -Bi(O) we 
have 

Aaw{x) = (A„)m(x^-""'') - {Aa)u{x) = - - F{u{x))] < 0, 

and we infer from Lemma 4.4 that {d / dxi)w{x) < for all x G Tx^^^ fl i?i(0). In 
conclusion, Amax < implies Amax £ A which by Step 2 forces sup A > Amax- This is 
a contradiction. 

4.3 Remarks and open questions 

The approach developed in the last sections invites to try to use the moving planes 
method also for solutions u of (3.1), that is to show radial symmetry about some 
point for the corresponding problem on all of R'^. In fact, this has been done in 
the case of the classical Laplacian starting with the work of Serrin [21]; see also 
[3, 9, 12, 14]. Observe that steps 2 and 3 from Section 4.2 can be carried over easily 
to the infinite space setting. In words, once we have, for any given normal direction, 
a "good" hyperplane T^, we can push it along until we hit a symmetry center. Here 
"good" means that A lies in the obvious analogue of (4.3). Unfortunately we have 
not been able to show that there are any good hyperplanes to start at all. Finding 
such a hyperplane would amount to proving that for A sufficiently negative we have 
u{x^) — u{x) > on T,\. This problem is of course trivial for u on a finite ball (see 
step 1 in Section 4.2), and it can be solved in the classical case by considering the 
differential inequality and boundary values one obtains for w{x) = u{x^) —u{x) on Ea. 
However due to the non-locality of Aq the question becomes more intricate in our case: 
Consider without loss of generality a hyperplane Tq through the origin and a w that 
has a reflection anti-symmetry w{x) = —w{x^) and solves Aaw{x) + ■ip{x)w{x) = 0, 
where ip is a positive function whose size and decay behavior can be controlled, given 
reasonable decay assumptions on the original u. If we could show that then w does 
not change sign in Eq, we would have found a "good" hyperplane. 

It is also tempting to conjecture that Theorem 4.1 holds if Aq, is replaced by the 
generator of a Levy process with rotationally symmetric increment distribution. 
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